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Abstract 

We simulate the formation and the evolution of global strings taking into account 
the expansion of the universe and the concomitant change of the effective poten¬ 
tial, that is, the change from the restoration stage of the global f/(l)-symmetry 
to the broken stage. Starting from the thermal equilibrium state, we run simula¬ 
tions for 100 sets of different initial conditions in 2-|-l dimensions. As a result, we 
find deviation from scale invariance. This is because the energy loss mechanism is 
pair annihilation rather than intercommutation, so that a string and an anti-string 
attract each other under a logarithmic potential. 


1 Introduction 


Since the possibility of the restoration of broken symmetries in the very early universe 
was suggested |l[], the cosmological consequences of the accompanying phase transitions 
have been investigated by a number of authors p. As Kibble hrst pointed out, topolog¬ 
ical defects of various types may be formed, depending on the pattern of the symmetry 
breaking (see also Ref. 0 for review). If the symmetry group G breaks into a subgroup 
H and the coset space G/H is not simply-connected, strings are formed. Strings as local 
defects have been extensively investigated in the context of structure formation theory 
On the other hand, global strings are predicted in the axion theory of the strong 
OF problem [^], and they may have radiated significant amount of axions as the universe 
cools down from the Peccei-Quinn breaking scale to the energy scale of QCD |]^. Also, in 
condensed-matter systems, analogous vortices are observed p|]. 

Most studies of the formation of strings have employed Monte-Carlo simulations. 
Vachaspati and Vilenkin hrst examined the initial conhguration of a string network. 
They found that it has the statistical properties of a Brownian random walk (see also 
Ref. [^), and found that a string network is made of closed loops with a scale-invariant 
size distribution and inhnite strings occupying 80% of the total string length. Since that 
time, many authors have reexamined the problem in more realistic situations than the 
original, focusing especially on the existence and the fraction (if any) of inhnite strings 
m- However, as seen from the fact that the above mentioned studies were carried out 
exclusively by the Monte-Carlo simulations, little work on the dynamical aspects of string 
formation has been done. 

One exception is Borrill’s paper 


12|, where the simulation of the bubble nucleation 


process in a hrst order phase transition is performed using a lattice-free method. However, 
in his approach the subsequent evolution of the string network could not be followed. In 
Ref. the Langevin-equation approach is considered. But in that study the ehect of 


the expansion of the universe was considered only partially. This is quiet diherent from 
the real dynamics. Ye and Brandenberger 0 , on the other hand, performed a numerical 
simulation properly taking cosmic expansion into account. But there, a Mexican-hat 
potential was adopted from the outset, and the simulation domain was so small that the 
whole system fell within one horizon before they could conclude that the network had 
reached the scaling regime. The evolution of local (gauged) strings has been investigated 
by using the Nambu-Goto action. There are many analytical and numerical [|l^ 
approaches, and it has been confirmed that the string networks relax to the scaling solution 
by intercommutation which generates closed loops that eventually disappear by radiating 
gravitational waves or particles. On the other hand, as for global strings, the Kalb- 
Ramond action must be used instead of the Nambu-goto action in order to examine 
the evolution. But in that approach it is difficult to include the long-range force between 
inhnite strings. Here, we use an approach based on a complex scalar held. 

In this paper, we simulate both the formation and evolution of the global string 
network by solving the dynamics of the complex scalar held directly. Starting from a 
symmetry-restored stage, we follow the dynamics of the complex scalar held, fully consid¬ 
ering the expansion of the universe and the associated change of the ehective potential, 
which exhibits spontaneous symmetry breaking below the critical temperature. 

The simulation is made in 2-1-1 dimensions as a hrst step to analyze the full 3-1-1- 
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dimensional analysis, because the spatial resolution and total cosmic time are limited due 
to our limited computing power. Note that in 2+1 dimensions strings are point-like, and 
each such string corresponds to an inhnite string in 3+1 dimensions so that the energy 
loss mechanism is pair-annihilation rather than loop-formation through intercommutation. 
The direct emission of goldstone modes from strings is suppressed because of the symmetry 
of the system. In a forth-coming paper we will report the results of our simulations in 


3+1 dimensions and their cosmological implications |18 


The remainder of this paper is organized as follows. In the next section, we establish 
the formulation for the numerical simulation. In §^, we give the result of the numerical 
simulation and judge whether during the evolution the scaling behavior is observed as 
expected from the analogy of local strings. In §^, the results obtained from the simulation 
are explained analytically using a phenomenological model. Finally, is devoted to 
summary and discussion. 


2 Formulation 


We consider the following Lagrangian density for a complex scalar held <h(x): 

£[$] = - i/[<F]. (1) 

Here, the potential H[<I>] is given by 

r[4>] = ^ (OK* - y)". ( 2 ) 

Then the one-loop hnite temperature effective potential in the high temperature 

limit is written as 

Hefr[$] = • (3) 

For T > Tc = y/^T], the potential Wfr has a minimum at <I> = 0, and the U{1) symmetry 
is restored. On the other hand, new minima |‘F|mm = V\/^ ~ (T/T^)'^ appear and the 
symmetry is broken for T < (Fig. |^). In this case the phase transition is of second 
order. 

In the expanding universe the effective Lagrangian density /Iefr[‘h] is given by 

- Wff[<I>] , (4) 


where is identihed with the Robertson-Walker metric. Then the equation of motion 

is 

4(x) + = -KbW , (5) 

a[tY 

where the prime represents the derivative d/d^\ and a{t) is a scale factor. The Hubble 
parameter H = d{t)/a{t) and the cosmic time t are given by 




Svr TT^ 

3M2i ^ 


9*T\ 


t 


1 




2^2 ’ 


( 6 ) 
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where Mpi is the Plank mass, is the total number of degrees of freedom for the relativistic 
particles, and radiation domination is assumed. We dehne the dimensionless parameter ( 
as 

T] yWTr^g^r]'^ J 


In our simulation, we take (^ = 10, which corresponds to rj (10^^ — 10^®) GeV, but the 
essential result is independent of this choice. The energy density of each point is written 
as 

pix) = + Kfr[^>] • (8) 

a{ty 

We simulate the system in 2+1 dimensions from the initial time ti = fc/4, that is, 
Tj = 2Tc, to the hnal time tf = 150+ = 37.5+, when Tj ~ T^/O.!. Since the G(l) 
symmetry is restored at the initial time t = ti, we adopt as the initial condition the 
thermal equilibrium state with the mass squared. 


m 


2 


d^V,s[m 

d\^\^ 


l-K+O 


(9) 


which is the inverse curvature of the potential at the origin at t = tj. In the thermal 
equilibrium state, $ and <h are Gaussian distributed with the correlation functions. 





The functions <I)(a;) and ^{y) are uncorrelated for x ^ y. We generate these helds for 
the initial condition in the momentum space, because the corresponding helds ^{k) and 
^{k) are uncorrelated there. 

Hereafter we measure the scalar held in units of t and x in units of ti, and the 
energy density in units of t~^. The equation of motion and the total energy density are 
given by 


4>(l) + - iv=<l.(a:) = - J>t _ (12) 

1 1 / 25 \^ 25 

p{x) = d>(a;)$t(a;) + _v<h(x) ■ V$t(a;) + _ (^|$|2 - -j + -|<I>P , (13) 

where A is set to unity.|^ The scale factor a{t) is normalized as a(l) = 1. We arrange 
512 X 512 lattice points for our simulations. The lattice spacing 6x is set to 0.1, and its 

^Strictly speaking, the high temperature approximation is valid only when A is much smaller than 
unity. However, a change of A only scales the physical quantity and is not relevant to the final result. 
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physical length fephys grows as a(t)6x = 0.1The time step 6t is chosen to be 0.01. 
At sufficiently low temperatures, the typical width of the string d is nearly 1/r]. Then, 

(~ 245 at t = tf) (14) 


(~ 0.98 at t = tf) (15) 

Even at the hnal time tf, the resolution of the simulation is sufficiently good, and the 
simulation box is larger than the horizon scale, so that we can safely impose a periodic 
boundary condition. 

Under these circumstances, we simulated the system using the second order leap-frog 
method and the Crank-Nicholson scheme over 100 runs with a different random realization 
of the initial conditions. 


H- 


Sx 


= 20 ^^/^ 


phys 


dXr- 




3 Results 


First we depict in Figs. how strings are generated as the symmetry is spontaneously 
broken and the gradient and kinetic energy density decreases. In these hgures, the total 
energy density at each lattice point is depicted at t = 1, 50,100,150. One can see that 
initially the held huctuates very strongly and the gradient and kinetic energy dominates 
over the potential energy. But as the universe expands and the temperature decreases, 
the false vacuum energy becomes dominant, and the string conhgurations become stable. 

Next, in order to judge whether the string network really enters a scaling regime in 
this case, we count the number of strings per horizon volume at various times. Since 
spacetime is discretized in our simulations, a point at which 4) = 0 corresponding to a 
string core is not necessarily situated at a lattice point. In the worst case, a point at 
which 4) = 0 lies at the center of a plaquette. Therefore in order to identify a string, we 
stipulate that a lattice is identihed with a part of a string if the potential energy density 
there is larger than that corresponding to the held value of a static string solution at 
r = 5xphys/ \/2,Q which is obtained by solving the equation 


d'^p Idp p 


VM 


0 , 


( 16 ) 


with 4)(r, 9) = p{r)e^^. Then we regard a connected region of lattices satisfying the above 
condition as one string. Around f = 20, the kinetic and the gradient energy is sufficiently 
suppressed in comparison with the vacuum energy of the string so that we can identify the 
string unambiguously. We carried out this identihcation procedure every 10 steps after 
t = 20. 

Figure ^ displays that time evolution of the number of strings per horizon volume 
averaged over 100 realizations. It is clear that the number is not constant but gradually 
increases with the cosmic time. Since the rate of increase is a decreasing function of 
time, one may think that we have observed only a relaxation stage and that it will 

^In fact, the potential energy is peaked so highly at the string that the identification of the string 
depends very little on this criterion. 
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eventually approach a constant. In order to judge if this is the case, we run a different 
set of simulations with different initial conditions, in which strings and anti-strings are 
arranged by turns every two lattices, like a checker-board, and <I> is set to 0 for all lattices. 
The initial string number per horizon volume in this case is as large as 100. The result 
is depicted in Fig. ^ together with the previous result. The number per horizon volume 
rapidly decreases once and then rebounds and continues to increase, as in the previous 
case. This result implies that the upward trend in the previous simulation is genuine, 
because, if it were spurious, it would gradually decrease and approach the asymptotic 
constant value from above. In the next section we present a phenomenological analytic 
model to explain why a departure from the scaling solution is obtained. 


4 Analytic explanation 


In the present simulations in 2-1-1 dimensions, strings are represented by point-like quasi¬ 
particles and anti-particles depending on the direction of the winding. Here, strings cannot 
intercommute with each other to form loops, but rather they can only pair-annihilate upon 
collision. Between a string and an anti-string works an attractive force proportional to 
the inverse separation, which is a consequence of the extended held gradient as a global 
defect. Seen from the simulation result, as the system relaxes sufficiently and enters a 
quasi-equilibrium state, a string usually remains at rest because it loses kinetic energy due 
to cosmic expansion. Once within the horizon, a pair of strings feel the above long-range 
force, and they start moving toward each other. Therefore, we have only to elucidate the 
dynamics of a pair of particles obeying such a force in an expanding universe. 

Let n{t) be the average number of strings per unit physical volume at the cosmic time 
t. When the system is in a quasi-equilibrium state, the probability per unit time P{t) 
that a string annihilates with an anti-string is written by 

Pit) ~ ^ . (17) 

where T{t) is the period required for a pair of strings at rest with mean separation, 
~ l/n{t), to pair-annihilate. Since n{t) is diluted by cosmic expansion and the above 
annihilation process, it satishes the Boltzmann equation. 


dn{t) 

dt 


—Pit) nit) — 2Hnit) , 
nit) nit) 


(18) 


In order to determine the time dependence of T(f), we investigate the relative motion 
of particles subject to an attractive force proportional to the inverse separation in an 
expanding universe, which is described by the equation 


d^r 1 dr 2 ti 

dp ^ t dt ait)‘^r ’ 


(19) 


where r is the comoving separation. 
it/ti)^^‘^r, this equation reads 


d'^x 


dP 


In terms of the physical separation x = a(f)r = 

lx 2 ti 

AP X 
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For the moment, we neglect the dissipation term due to cosmic expansion. This is justihed 
below. Then integration of Eq. (pO]) over t gives 


1 / dx\ 

2 \dt j 


+ 27 rln f—) 
Vxo/ 


0 , 


( 21 ) 


where % = 0 and x = Xq ai t = 0. Then the period T is estimated as 



where d is the core radius. Before inserting this into Eq. ([T8|) to solve for n{t), we justify 
the omission of the dissipation term. A pair of strings at rest initially with separation xq 
approach a separation x in a time interval t given by 


t ~ 



(23) 


Then the dissipation term, x/(4t^), reads 7rx/(4xo ln(xo/x)). Since x is smaller than xq, 
the term from the potential force in Eq. (^0|), 27r/x, always dominates over the dissipation 
term. 

Taking xq as the mean separation, 1/ \Jn(t), Eq. (^) becomes 


dn(t) 

dt 



(24) 


where l/yn{t) is nearly proportional to the cosmic time from the simulation result, and 
the proportional coefficient is A. Substituting f{t) for n{t)t, we have 


1 df _ 1 



Integrating over both sides gives 



(25) 


( 26 ) 
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where B represents a numerical factor ~ 0{1) and C is an integration constant determined 
from the initial conditions. 

Then n{t) is given by 


n{t) ~ 


D 


1 + 


■Tt) 


(27) 


where D = ~ d/A ~ d/\/lD. Therefore the number of strings per horizon 

volume is given by 


n(t)H ^ 


AD 


>■'(t) 


1 + 


cr>i/2 

P/2 


In 


n 2 


AD In 


(asymptotically) 


(28) 


Taking D 0.1, E rsj 0.4 (^d 1.0) and C —0.09, we can ht the results of the simulation 
by the above formula. This is depicted in Fig. |^. The number of strings per horizon volume 
is not a constant, but proportional to (Int) asymptotically. Note that only the constant 
C depends on the initial conditions. Taking C ~ 0.8 with almost the same D and E, that 
from the checker-board initial conditions can also be htted. This implies that our htting 
formula is consistent. 


5 Summary and Discussion 

In this paper, we have reported the dynamics of a complex scalar held from the symmetric 
phase to the broken phase in 2-1-1 dimensions. Starting from the thermal equilibrium state, 
we have observed how strings are formed as the temperature decreases and the symmetry 
is spontaneously broken. We also hnd that, as the string system grows, it does not go 
into a scaling regime, where the number of strings per horizon volume stays constant, but 
the number increases in proportion to (Int) asymptotically. This can be interpreted as 
follows. Our simulation was performed in 2-1-1 dimensions and strings are in fact point¬ 
like. Thus, the dominant energy loss process is not the formation of loops through the 
intercommutation of inhnite strings, but rather pair-annihilation of straight strings. They 
move under logarithmic potentials, which causes deviation from scale-invariance. 
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Veff[0] T>Tc T=Tc T = 0 



Figure 1; One-loop finite temperature effective potential of the complex scalar 

field. 
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Figure 2: The total energy density for each lattice point is depicted aX t = U = 1. 
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Figure 3: That at t = 50. 
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Figure 4: That at t = 100. 


13 

















time 


Figure 6: The number of strings per horizon volume is depicted. Dots represent the 
average value over 100 realizations. Error-bars the standard deviation. 
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Figure 7: Open dots represent the number of strings per horizon volume averaged over 
those from 100 initial conditions which are the thermal equilibrium states. On the other 
hand, hlled dots from checker-board initial conditions. Also, a solid line represents the 
htting formula for simulation results from thermal initial conditions and a dotted line 
from checker-board ones. 
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